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1. Introduction 

This note is on inverse spectral theory for the Schrodinger operator on a flat 
two-dimensional torus with electric and magnetic potentials. This problem can 
be remarkably rigid. For generic flat tori, if the variation of the magnetic field is 
strictly less than its mean, and the total magnetic flux on the torus is ±27r, then the 
spectrum of the Schrodinger operator determines both the electric and magnetic 
fields. This is in marked contrast to both the Schrodinger operator without a 
magnetic field (see [3]) and the case of a magnetic field of mean zero (see [1]). 
In both those problems there are large families of isospectral fields, and rigidity 
results are much more difficult to obtain (see also [2] ) . The observation that there 
can be spectral rigidity when the total flux is ±27r is due to Guillemin ([5]). Here we 
give a short proof of the slightly stronger result stated above. Instead of thinking 
of the Hamiltonian as acting on functions with values in a line bundle over the 
torus RV-^: 

we think of the Hamiltonian as acting on functions on which are 
invariant with respect to the "magnetic translations" associated to L. However, 
these two settings are completely equivalent. Our assumption that the variation 
of the magnetic field B{x) is strictly less than its mean bo takes the simple form 
\B{x) - 6o| < l^ol for all x. 

The spectrum of the Laplacian plus lower order perturbations on flat tori 
has the feature that there are large families of spectral invariants corresponding 
to sets of geodesies with a fixed length. In analogy with results on S'^ Guillemin 
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proposed the name "band invariants" for these famiUes. The nice feature of the 
problem discusscid here is that only the simplest of the band invariants are needed 
to prove rigidity. 

The first complete solution of an inverse spectral problem was Mark Krein's 
definitive analysis of the "weighted string", [9], [10]. Since that time many other 
inverse spectral problems in one space dimension have been solved (see [11]). In 
higher dimensions it is widely believed that, modulo natural symmetries and de- 
formations like gauge transformation, most problems will be spectrally rigid. How- 
ever, so far there have been relatively few settings where this has been proven (for 
instance those in [6] and [14]) and many interesting examples where it fails (see [12] 
and [4]). This should remain an active field of research for many years to come, and 
one can reasonably say that it began with the work of Mark Grigor'evich Krein. 

2. Proof of Gtiillemin's Theorem 

We begin with the smooth magnetic field B, periodic with respect to the lattice 
L in two dimensions, expanded in a Fourier series in terms of the dual lattice L* 

B{x) = 6/3e2^'^-^. 

For this magnetic field we introduce the magnetic potential A = {Ai,A2) with 
dx2Ai — 9x1^2 = B, chosen to be as periodic as possible, i.e. 

A = A° + A' = ^{X2, -xi) + b0e^-'^-^p2, -Pi)27ri{pf + ^|). 

/3eL»\o 

We also have a mean zero periodic electric field which is the gradient of the mean 
zero periodic potential 

V{x)= Y ^'/3e^"'''"'- 

The quantum Hamiltonian for an electron in these fields (with all physical con- 
stants set to 1) is 

H = {idx + Af + V. 

Let _D be a fundamental domain for L. To define the domain of H as an 
operator in L'^{D) we will use "magnetic translation operators" (see [13]). Letting 
{61,62} and {6^,62} be a basis for L and the corresponding dual basis for L*, 
define for linearly independent vectors Vi and V2 

Tju{x) = 6'''^-^u(a; + ej), j = 1, 2. 

Then the commutator [Ti,T2] is given by 

[n,T2]u{x) = (e*''^-^^ -6'''i-'==)e^(''i+''^)-^w(a; + 6i +62), 
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and the periodicity of and V implies that the commutator [H, Tj] is given by 

[H,Tj]u{x) = e*"^-"((z9, + A{x) + A"{e,))^ - {id., + A{x) ~ Vjf)u{x + e^). 

Thus, in order for the Tj's to commute with H we require Vj — — ^"(ej), and in 
order for the Tj's to commute with each other we require A^(ei)-e2 = —A^{e2)-ei = 
111 for some integer Note that this implies A^{e\) = 77/62 and AP{e2) = —nlel 
and 27r|/| = |&o|Area (D), and &oArea {D) = B{x)dx is the total magnetic flux. 
Hence the assumption &o 7^ is equivalent to nonzero flux, and it implies / ^ 0. 
Defining the domain of H to be the subspace of H^{R^) such that Tju = u, j = 
1,2, we make H a self- adjoint operator in L'^{D). 

As in many previous works we will look for spectral invariants for H by study- 
ing the wave trace. Letting y, t) be the distribution kernel for the fundamental 
solution for the initial value problem 

utt + Hu = QinRlx Mt, M(a;, 0) = /(x), 0) = 0, 

the distribution kernel for the corresponding initial value problem \a D xM.t is 

ED[x,v,t)^ TrT^E{x,y,t), (2.1) 

(m,n)GZ2 

where the operators Tj act on the x variable. Note that, since the principal part of 
+H is — A, E{x, y,t) =0 when \x — y]"^ > and the sum in (|2.1|) is has only 
a finite number of nonzero terms for i in a bounded interval. Thus \Ti,T2] — 
implies TjED{x,y,t) = Eu{x,y,t), j = 1,2. The fundamental spectral invariant 
for this problem is the distribution trace of the operator Ejj (i) corresponding to the 
kernel E]j{x,y,t). Conventionally (with all terms to be interpreted in distribution 
sense) this is written 

Tr{t) = I ED{x,x,t)dx. 
Jd 

To avoid degeneracies in the contributions to Tr{t) from the terms in (1), we 
assume that vectors in L have distinct lengths, i.e. 

d,d' € L and \d\ = \d'\ implies d = ±d' . 

Since E{x,y,t) is singular as a distribution in {x,y) only when \x — j/p = i^, it 
now follows that the singularity of Tr(t) at t — \mei + ne2\ comes from just two 
terms 

f [Tl^T^Eix, X, t) + Tf ™T2~"£;(a;, x, t)]dx. (2.2) 
Jd 

To determine the spectral invariants coming from the leading terms in the expan- 
sion of this singularity it is convenient to use the Hadamard-Hormander expansion 
[7], [8] for i?(x, y, t). Beginning with the forward fundamental solution, E^, defined 
by {dj + H)E+ = 5{t)5{x - y) and E+^{) for t < one has 

00 

E+{t, x,y) ^Y a^{x, y)e^{t, \x - y\) (2.3) 

i/=0 
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where Ci, is chosen so that [d^ — A)eiy = i^e^_i for v > Q and eo(t, \x — y\) is the 
forward fundamental solution for df — A. In two space dimensions this means 

e.(i, \x - v\) = 2--'^-'n-^l-'X'i-"\e -\x- yp) 

for i > 0, = for t < 0. For a > — 1 the distribution A"" is defined by 
X^{s) = (r(a + l))-is'' for s > and A'^(s) = for s < 0. Hence the coefficients 
an are determined by the recursion 

uanU + {x - y) ■ dxay - iA{x) ■ {x - y)ai, + Ha^^i = 0, 

where H acts in the variable x. Solving this with the requirement that ao(2/, y) = 1, 
we have 

ao(x, y) = exp(i / {x - y) ■ A{y + s{x - ?;))fis) and (2.4) 
Jo 

ai{x,y) ^ -ao{x,y){j V{y + s{x ~ y))ds + b{x,y), (2.5) 
Jo 

where b{x,y) is determined by A{x). The fundamental solution E{x,y,t) is given 

by 

i;(a;, y, t) = at(S+(i, x, y) - S+(-t, x, y)). (2.6) 

We define 

I[d)= [ f,-'-A°{d)-x+if^] d-Aix+sd)ds^^ ^ r ^i(2A"(x)-d-J^d-A\x+sd)ds)^^^ 
J D J D 

and 

J{d)= I [f {V{x + sd) + b{x + sd,x))dsy^^^°^'''>-'^~^od-AHx+sd)ds)^^^ 
Jd Jo 

From (|2.2p - (|2.6p one sees that I{d)+I{—d) and J{d) + J{—d) are spectral invariants 
for H. However, the periodicity implies that I{d) — I{—d) and J{d) — J{—d). 

The rest of this article is devoted to studying I{d) and J(d). We have d = 
mei + ne2 = k{moei + 71062), fc S N and gcd(mo, no)=l. Let S = —noel + 171062- 
Then we have 

A°(d) = ^(^d2,-di)=7rklS. 

Since e'^^'^^^''^ds = when (3-d ^ 0, the terms in the Fourier series for A^ which 
contribute to 1(d) have (3 - d = 0, and this implies 

/3-P'5 = |^(d2,-di), pez\0. 

Hence, d ■ {(32, -(3i){2ni{(3l + ^|))~^ = ikl{pba)~^ , and I{d) reduces to 
/ exp(27rzfcZ(-<5-a;+ V ^^e'^'^'-^))dx. 



Remark on spectral rigidity for magnetic Schrodinger operators. 



5 



Defining 

pez\o pGZ\o 
(note that j^Al{s) = Bs{s)), we have 

I{d) = / exp{~i2TTkl{S ■ x + ^As{S ■ x))dx. 
Jd bo 

Extending 5 to a basis for L*, {6,6'}, and letting {7,7'} be the dual basis for L, 
we make the change of variables x — sj + uj' , and choose 

D = {s-f + u-f' : < s,u < 1}. 

Then we have 

/■I 1 

I{d) = c{d) / exp{~2TTikl{s + —A]{s))ds, 
Jo bo 

where c{d) is the Jacobian factor, and only depends on d. Since we have this 
spectral invariant for all fc 7^ 0, it follows that 

f' f{s + lAl{s))ds (2.7) 
Jo "0 

is a spectral invariant for any function / which can be expanded in terms of 
{e-2^''='!^}fcez, i.e. for any / e Lf^^{R) which has period l/l. 

Theorem 2.1. Assume that I — 1 and |6o| > max|_B(a;) — 6o|. Then the spectrum 
of H determines B. 

Remark 2.2. Since 60 is the average of B(x) on a fundamental domain, the hy- 
pothesis here is a constraint on how much B varies instead of constraint on the 
size of B. 

Remark 2.3. Since we assume 1 = / = A°(ei) • e2/7r — |6o|Area(Z))/27r for a 
fundamental domain D, this assumption fixes |6o| when L is fixed. 

Proof of Theorem \2.1\ Since B^{x) — ^^{B{x + sd) — bo)ds, the hypotheses imply 
that the derivative of s + Ag{s)/bo is strictly positive and the inverse function s{y) 
to y = s + Al{s)/bo is defined on the range of s + Aj(s)/&o for s G [0, 1]. Since A] 
has period 1, the range is / = [Aj(0)/6o, Al{0)/bo + 1]. Letting / in (|2T7l) tend to 
the 5-function at y, the limit of ()2.7p is s'{y) ii y = y{s) for s G [0, 1]. If y ^ y{s) 
for s G [0, 1] then the hmit of (|2.7p is s'(y*), where y* G /, and y* = y mod 1. 
In other words taking these limits we recover a function of period 1 in y which 
agrees with s'{y) on /. Thus we recover Ag(s) modulo an additive constant, and 
we obtain Bs{s) by taking the derivative. Since we can carry out this argument 
for all prime elements 6 £ L*, we recover the full Fourier expansion of B. □ 
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We now turn to the recovery of V. The preceding analysis shows that, keeping 
the same d G i as above, the spectral invariant J{d), modulo terms determined 
by A{x), reduces to 

j{d)^c{d)[ l^5(s)e-2"^'=(^+^^(^»ds, (2.8) 
Jo 

where 

pez\o 

This immediately gives the following: 

Theorem 2.4. Under the hypotheses of Theorem \2.1\ the spectrum of H determines 
V. 

Proof of Theorem \2.4\ Since we are assuming the hypotheses of Theorem 12.11 we 
have the function s{y) and can make the substitution s = s{y) in (j2.8l) . That gives 

J{d) = c{d) / F5(s(y)e-2-fc!^s'(y)dy, 

J Alio) 

but, since y{s + 1) = y{s), we can extend s{y) smoothly to the whole line by 
defining s(y + 1) = s{y) + 1. Thus, since Vs{s) has period 1 in s, we have 

J{d)^c{d) fvs{s{y))s'{y)e'^^'^^ydy. 
Jo 

Since we have this spectral invariant for k G Z\0, we recover the Fourier series 
of Vs{s{y))s'{y), and, hence, since s{y) is determined by Al{s), we have Vs{s). As 
before, since we can carry out this argument for all prime elements S d L*, we 
recover the full Fourier expansion of V^. □ 

Remark 2.5. If / ~ p/q, p,q E N, for the lattice L, then I — 1 for the lattice Lq 
generated by ciei + 0262 and dici + d2 when p{cid2 — 02^1) = q. So if B{x) and 
V{x) are periodic with respect to Lq, Theorems 1 and 2 apply in the sense that 
the spectrum of H on the torus M^/£o determines B{x) and V{x). Note that B{x) 
and V{x) will automatically be periodic with respect to Lq when / = 1/q. 
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